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Abstract: The modelling of human exposure in realistic exposure scenarios is complex, because 
several parameters (e.g., the source design, the frequency band, the orientation of incident fields, 
the morphology and posture of subjects) vary and influence the dose. Deterministic dosimetry, so 
far used to quantify human exposure to electromagnetic fields (EMF), is highly time consuming 1f 
the variations of those parameters are considered. Stochastic dosimetry is an alternative approach to 
assess EMF exposure and consists in building analytical approximations of the exposure at a 
parsimonious computational cost. In this study, it was used to assess the influence of magnetic flux 
density (B) orientation on fetal exposure at 50 Hz using the polynomial chaos (PC) theory. A PC 
expansion of induced electric field (E) in each fetal tissue at 7 months of gestational age (GA) was 
built as a function of B orientation. Maximum E in each fetal tissue was estimated for different 
exposure configurations and compared with the limits of the International Commission of Non- 
Ionising Radiation Protection (ICNIRP) Guidelines 2010. PC theory resulted in an efficient method 
to build accurate approximations of E in each fetal tissue. B orientation influenced E with a 
variability across tissues in the range from 10% to 25% with respect to the mean value. However, 
varying B orientation, maximum E in each fetal tissue was below the limits of ICNIRP 2010. 
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INTRODUCTION: VARIABILITY OF HUMAN EXPOSURE TO EMF 


In the evaluation of the human exposure to EMF it is necessary to take into 
account several parameters, that vary in a real exposure scenario: 


> Source (location, design, frequency) 
> Orientation of the incident field 

> Environment 

> Morphology 

> Dielectric properties 


> Posture 
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INTRODUCTION: STUDY OF EMF EXPOSURE 


DETERMINISTIC DOSIMETRY STOCHASTIC DOSIMETRY 





Orienta ti Induced electric fields Orientation i pelli 
Morphology Pietre orren tdeneily ROSI SLOEY e Electric current density 
Posture o SAR e Posture è SAR 

° Dielectric properties 


Dielectric properties 





HIGHLY TIME CONSUMING! 


Polynomial Chaos (PC) decomposition: 
approximation of Y on a suitable basis of 
orthogonal polynomials (Wiener 1938; Xiu et 
al., 2002) 
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INTRODUCTION: POLYNOMIAL CHAOS (PC) PRINCIPLE (1) 
Soize and Ghanem, 2004 


Hp: E[Y?]<+c0 
polynomial 





deterministic coefficient 


Construction of the polynomial basis W(X) 


X= vector of K indipendent input parameters each characterized by probability density 
function (PDF) fy; 


K 
‘> Wa (X) — BEZ — Ha, (X,) XK wees Xa, (Xx) polynomial 
J=1 
Ia; = family of polynomials orthogonal respect to each f,, 


Distribution Support Polynomial 
a;= maximum degree of la; -a a a aa 


Gaussian R Hermite 
, Uniform [—1,1] Legendre 
polynomial degree: |a|=a,+...+a,< p 
Gamma (0,+00) Laguerre 
W(X) size: P= (K+p)!/ K!p! Chebyshev (—1,1) Chebychev 
Beta (—1,1) Jacobi 
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INTRODUCTION: POLYNOMIAL CHAOS (PC) PRINCIPLE (2) 


Experimental Design 
K={x'2) (2) X(N) (N+) 


Estimation of a Observations l 
y=fy(1) y(2),.. yl) y(n } 








Y = X ai pi(X) 


i=0 


NO 


Validation set 
Vial 
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INTRODUCTION: POLYNOMIAL CHAOS (PC) PRINCIPLE (3) 








Experimental Design 
M=(xl2}, 2612), ...,x(N1} 





Least Angle Regression algorithm (LAR, Efhron r= Y aoo 
et al., 2004) adapted by Blatman et al., 2011 to 
the PC theory 


Estimation of coefficients of PC 
expansions and collection of PC 
expansions 


Validation set 
Yyal 





1. LAR generates a collection of PC expansions (the first expansion includes a single 
polynomial ,, the second two polynomials (W,, ,), until m=min(P,N-1)) 


2. LAR chooses the best PC expansion by Leave-one-out cross-validation. 
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OBJECTIVES 





Estimation of induced E field in each fetal tissue exposed to ELF-MF at 50 Hz 
changing the B-field orientation by means of PC decomposition: 


1. Build a PC expansion of E in each fetal tissue; 
2. Statistical analysis of fetal exposure; 


3. Identification of the worst-case exposure scenario with respect to the 
ICNIRP Guidelines 2010. 
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MATERIAL AND METHODS: PC APPLIED TO ELF-MF PROBLEM 





Y = M(X) = Xo aj W(X) 


e Y=99th percentile of induced E in each fetal tissue 


6 uniform distribution [0,180° ] 
e X input random vector: 
o uniform distribution [-180° ,180° 4 


& 


¢ W(X) polynomial basis: Legendre polynomials 


Distribution Support Polynomial 


Gaussian R Hermite 
Uniform [-1, 1] Legendre 
amma D, +00 aguerre 


Chebyshev (—1,1) Chebychev 
Beta (-1,1) Jacobi Xiu et al., 2002 
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MATERIAL AND METHODS: EXPERIMENTAL DESIGN 






Poi 
Y= >: ajwp;(X) 
j=0 


Estimation of a. 





e Pregnant woman model at 7 months GA 
(provided by IT’IS Foundation) 


Validation set 
Yvat 


e Low Frequency Solver SEMCAD X 





e Exposure: 200 uT uniform MF at 50 Hz 


e Simulation setting adopted in Liorni et al., 2014 





Experimental Design: generation of N couples (6; «) by Quasi-Monte Carlo 
method based on Sobol’s function depending on joint PDF fy: 
b, Din 


b 
Du DLL L 
ues 2 tzt t Smr 


Observations: estimation of N induced E in each fetal tissue on the 
experimental design. 
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MATERIAL AND METHODS: VALIDATION OF THE PC EXPANSION 


Experimental Design 
K={xf2} x !2), al} 


Deterministic dosimetry 








Observations 
y={y! 1) yl? l, aa y} 


PROCEDURE OF VALIDATION 

1. S couples of O and ọ different from the experimental design X 
2. Estimation of E on S by deterministic dosimetry (y,„„P) 

3. Estimation of E on S by PC expansion (y,,)°°) 


4. Calculation of the percentage mean square error (pMSE) between y,’ 
and y,,\°° 
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RESULTS: CHOICE OF PC EXPANSION 





Example: Trend of pMSE increasing the number of observations to build a PC expansion of 
E induced in the fetus whole-body. 


1.8 1 
1.6 | 
1.4 | 
Lez 

TE 
0.8 7 


pMSE (%) 


0.6 | 
0.4 | 
0.2 47 





0 
50 100 150 200 250 300 


Number of Observations 


The validation has been performed for each fetal tissue at 7 mGA. N=300 


observations is suitable to build all PC expansions with pMSE not higher 
than 0.17%. 
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RESULTS: STATISTICAL ANALYSIS OF FETAL EXPOSURE 


Estimation of mean and standard deviation from PC coefficients (Blatman et al., 2011) 








9 , 
7 MGA 
9 È 
7 $ 
6 $ 
æ 9 TT 
E 
= 
È 4 
us 
3 $ 
2 $ 
1 + 
yea 
OO È wk oe adh È - LC NAN è @& © CS ~ Oy 
PEOL O Po a EEK FT E FT È 
Tere Ye © SR go N 


e Fetal skin, fat, liver, SAT present the highest mean E (up to 5.09 mV/m); 
e Mean E induced in bone tissue up to 3.08 mV/m; 


e Variation of E higher than 10% and up to 25.3% in gallbladder. 
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RESULTS: WORST-CASE SCENARIO RESPECT TO ICNIRP 2010 





Ews E, is the max E among the 10000 random generated in each tissue 


% 





WS% = 
lim 


Em ICNIRP 2010 Basic Restriction for the General Public at 50 Hz 


e CNS of the head (limit of 0.02 V/m): 23% of the limit 


e All the other tissues of head and body (limit of 0.4 V/m): 


e Max E always under the limit, 
with WS% lower than 2%; 


e Highest max E found in the fetal 
tissues with also the highest 
mean E. 
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CONCLUSIONS 


Polynomial Chaos decomposition results an efficient method to study the variation of human 
exposure to EMF; 


e Study of the variation of B-field orientation at 50 Hz on fetal exposure: 


v The highest mean and maximum exposure found in fetal skin, fat, SAT, liver at 7 months 
gestational age; 


vY All the fetal tissues always result under the limit of the ICNIRP 2010 for the General 
Public at 50 Hz; 


v E distribution in some tissues is significantly influenced (up to 25%) by B variation. 
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